Josephson effect for SU(4) carbon nanotube quantum dots 
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We present the theory of the Josephson effect in nanotube dots where an SU (4) symmetry can be 
realized. We find a remarkably rich phase diagram that significantly differs from the SU{2) case. In 
particular, vr-junction behavior is largely suppressed. We analytically obtain the Josephson current 
in various parameter regions: (i) in the Kondo regime, covering the full crossover from 5'f/(4) to 
SU{2), (ii) for weak tunnel couplings, and (iii) for large BCS gap. The transition between these 
regions is studied numerically. 
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Several experimental groups have recently 
started to study the Josephson effect in ultra-small 
nanostructuresfi where the supercurrent can be tuned 
via the gate voltage dependence of the electronic levels 
of tlie nanostructure. An important system class 
where supercurrents have been successfully observed^ 
is provided by carbon nanotube (CNT) quantum dots. 
In many cases, the experimental results compare quite 
well to predictions based on modeling the CNT dot 
as a spin-degenerate electronic level with SU{2) spin 
symmetry, where the presence of a repulsive on-dot 
charging energy U may allow for a (normal-state) 
Kondo effect. Depending on the ratio Tk/^, where 
A is the energy gap in the superconducting electrodes 
and Tk the Kondo temperature, theory^"— predicts a 
transition between a unitary (maximum) Josephson 
current for A <g; Tr-, possible thanks to the survival 
of the Kondo resonance in that limit, and a 7r-junction 
regime for A 3> Tk, where the critical current is small 
and negative, i.e., the junction free energy F{Lp) has a 
minimum at phase difference = tt as opposed to the 
more common 0-j unction behavior. 

Recent progress has paved the way for the fabrica- 
tion of very clean CNTs, resulting in a new generation 
of quantum transport experiments and thereby revealing 
interesting physics, e.g., spin-orbit coupling effects^ or in- 
cipient Wigner crystal behavior In ultra-clean CNTs, 
the orbital degree of freedom {a = ±) reflecting clock- 
wise and anti-clockwise motion around the CNT circum- 
ference (i.e., the two K points) is approximately con- 
served when electrons enter or leave the doti^i Due to 
the combined presence of this orbital "pseudo-spin" (de- 
noted in the following by T) and the true electronic spin 
{S), an enlarged S'C/(4) symmetry group can be real- 
ized. In addition, a purely orbital SU{2) symmetry arises 
when a Zeeman field is applied. Experimental support 
for this scenario has already been pub lished^ (for the 
case of semiconductor dots, see Ref. [l^), and several as- 
pects have been addressed theoreticallyJii^ In particu- 
lar, the SU{4:) Kondo regime is characterized by an en- 
hanced Kondo temperature and exotic local Fermi liq- 
uid behavior, where the Kondo resonance is asymmetric 



with respect to the Fermi level. However, so far both 
experiment and theory have only studied the case of 
normal-conducting leads, where conventional linear re- 
sponse transport measurements cannot reliably distin- 
guish the SU{4:) from the SU (2) scenario Here we pro- 
vide the first theoretical study of the Josephson effect 
for interacting quantum dots with (approximate) SU (4) 
symmetry, and find drastic differences compared to the 
standard SU(2) picture. In the Kondo limit, a quali- 
tatively different current-phase relation (CPR) is found, 
with the critical current smaller by a factor « 0.59. The 
usual TT-j unction behavior is largely suppressed, but new 
phases do appear and time-reversal symmetry can be 
spontaneously broken. Our predictions can be tested us- 
ing state-of-the-art experimental setups, and offer clear 
signatures of the SU{4) symmetry in very clean CNT 
quantum dots. 

Model and formal solution. — We study a quantum 
dot (Hd) contacted via a standard tunneling Hamiltonian 
(Ht) to two identical superconducting electrodes {H^/^), 
H = Hd + Ht + Hl + Hfj. We assume that the dot has 
a spin- and orbital-degenerate electronic level e^cr = e 
with identical intra- and inter-orbital charging energy 
Hd = en + Un{n ~ l)/2 with fi = Eca4.rfo., 
where creates a dot electron with spin a 1= ± 
and orbital pseudo-spin projection a. Since the a = ± 
states are related by time-reversal symmetry (clockwise 
and anti-clockwise states are exchanged), we take the 
lead Hamiltonian as 
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creates an electron with wavevector fe 
- L/R, and is the single-particle 
energy. The tunneling Hamiltonian is Hf = 

J2jka,aa' (^'^aa' + ^X.-a') c]f^o.J^'a + ^^'C-, whcrC t (t) 

describes orbital (non-)conserving tunneling processes. 
Following standard steps,— the noninteracting lead 
fermions can now be integrated out. The partition func- 
tion Z{ip) = e~^^^^^ at inverse temperature (3 then reads 
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(we often set e = h = 1) 

Ziip) = Tr, (^e-^H,^e- ^0 drdr' DUrmr-r')D(r')^ ^ 

where the trace extends over the dot Hilbert space, T 
denotes time ordering, and we use the Nambu bispinor 
D ~ (c?g^, dg^, d^^, djj^) with even/odd hnear combina- 
tions of the orbital states, d^a — (rf+.o- + c?-,o-)/"\/2 
and doer = cr(d+,cr — d-^a)/V^- In this basis, the 
self-energy I](t) representing the BCS leads is diago- 
nal in orbital space. With the orbital mixing angle 
9 ~ 2tan~^(f/i) and the normal-state density of states 
1^0 = the even/odd channels are character- 

ized by the hybridization widths r^=e,o = (1 i sm6) T 
with r = TTi'o(t'^ + t"^). In what follows, we study the 
zcro-tcmperature limit and assume the wide-band limiti 
for the leads. The Fourier transformed self-energy is 
then expressed in terms of the 2x2 Nambu matrices 
r„ / —iuj A cos 7 
* -iuj 



A cos ^ 



The result 



P]) will now be examined in several limits. We start with 
the strong-correlation limit U oo, and later address 
the case of finite U. Note that Eq. ([IJ for 9 = corre- 
sponds to the SU{4:) symmetric case while for 9 = tt/2 
there is only one conducting channel with non-zero trans- 
mission which, under certain conditions, corresponds to 
the usual SU{2) model. 
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FIG. 1: (Color online) Josephson CPR in the Kondo limit 
for various Q. The Sf/(4) case corresponds to 61 = 0, the 
SU{2) case to 6 = n/2. The supercurrent is given in 
units of the unitary limit Ic = eA/h. 

Deep Kondo limit. — Let us first discuss the Kondo 
limit Tk S> a in the quarter-filled case, e < and 
{h) « 1. The Kondo temperature is given by Tk = 
£>exp(7re/4r)ii with bandwidth D. As in the SU{2) 
case^ the Josephson current at T = can be com- 
puted from local Fermi liquid theory, either using phase 
shift arguments or an equivalent mean-field slave-boson 
treatment)^ The latter approach yields the self-consistent 



dot level e and thereby the transmission probability for 
channel u = e,o^ 



(1± sin 6I)2t|. 
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In the SU{4) case (9 = 0), we have Te = To = 1/2, while 
the SU{2) limit {9 = 7r/2) has a decoupled odd channel, 
Te = I and To = 0- The CPR covering the crossover from 
the SU{4:) to the SU{2) Kondo regime then follows as 
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The known SU{2) result^ is recovered for 9 = tt/2. The 
5t/(4) CPR has a completely different shape, as shown in 
Fig. [TJ We note that the critical current = max[/((^)] 
is suppressed by the factor 2 — \/2 « 0.59 relative to 
the unitary limit eA/h reached for the SU{2) dot. The 
Josephson current in the deep Kondo regime is thus very 
sensitive to the 5J7(4) vs SU{2) symmetry. 

Perturbation theory in T . — Next we address the op- 
posite limit of very small T <|; A, where lowest-order 
perturbation theory in F applies. After some algebra, 
Eq. ([TJ for = yields the CPR of a tunnel junction, 
I{^p) = Ic sin((/3), where the critical current is 

/, = (4e(6)-e(-e))F(|6|/A)/o, (4) 

with the Heaviside function 8, the current scale /q = 
A(F/7rA)2, and (see also Ref.[ll) 



F{x) = 



{^/2f{l-x) 



arccos a; 
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In this L/ — > oo limit, the dot contains one electron for (fi- 
nite) e < 0, and thus wc have spin S = 1/2. Equation ^ 
shows that such a magnetic junction displays a 7r-phase. 
For the SU{4) case, the ratio /c(-|e|)//c(|e|) = -1/4 
is twice smaller than in the SU{2) case, i.e., 7r-junction 
behavior tends to be suppressed. This tendency is also 
confirmed for U A (see below), where the 7r-phase is 
in fact essentially absent. The factor 1/4 can be under- 
stood in simple terms by counting the number of possible 
processes leading to a Cooper-pair transfer through the 
doti ^^i^^ . When e > 0, there are four possibilities cor- 
responding to the quantum numbers (a, a) of the first 
electron entering the dot. However, for e < there is 
only one possibility since an electron already occupies 
the dot, and then only one specific choice of {a, a) allows 
for Cooper pair tunneling. This argument is readily gen- 
eralized to the SU{2N) case, where the above ratio of 
critical currents is obtained as — 1/2A^. 

Effective Hamiltonian for A ^ oo. — The partition 
function ([1]) simplifies considerably when A exceeds all 
other energy scales of interest. Then the dynamics is al- 
ways confined to the subgap region (Andrcev states), and 
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FIG. 2: (Color online) Phase diagram for A — !> oo. White 
regions correspond to {S, T) = 0, and green regions to 
{S,T) — 1/2. In the black regions, the ground state has 
{S, T) = for V? = and (S, T) = 1/2 for ip = n. For the 
blue region, we have (S, T) = at </3 = and (S, T) — 1 

quasiparticle tunneling processes from the leads (contin- 
uum states) are negligible. In particular, this allows to 
study the case U A. In fact, for A — )■ oo, with the 
Cooper pair operators b\ ~ dl^dl^ and ~ '^l-fC^e^i 
Eq. Ill) is equivalcntly described by the effective dot 
Hamiltonian 



(Similarly, the subspaces corresponding to Hj^gT)=i/2 
are obtained by letting dj,^ — > d\_^ and b\ — > b\.) 
With Tf. > To, the lowest energy is E^s.t)=i/2 = 

(^E^ +Es- [(i?3 - + ATI cos'{ip/2)] /2. (iii) 

Finally, the {S,T) — (1,0) sector is spanned by the 
two uncoupled two-particle states (i^g.(ij^|0), with 
(/3-independent energy Eg^i x^o = ^'2- In addition, 
there are two decoupled {S,T) = (0,1) states d^^dl^lO) 
with the same energy £'2- In the limit A ^ 00, this 
{S,T) = 1 sector is energetically unfavorable except 
possibly at = TT. 




-8/r 



H^ = Hd + cos(</j/2) [r,bi + rob2 + h.c] . (5) 

The resulting Hilbert space can be decomposed into 
three decoupled sectors^ according to spin S and or- 
bital pseudo-spin T (notice that these quantities are lo- 
calized on the dot for A 00). The ground-state en- 
ergy Eg{ip) — mm{E(^sT)) then determines the Joseph- 
son current I{ip) = 2d^Eg[Lp). (i) The {S,T) = sector 

is spanned by the four states {|0), 5| |0), 6j|0), 5|4|0)}, 
where |0) is the empty dot state. The matrix represen- 
tation reads 

/ Fecosf To cos f \ 
„ _ I Te cos f E2 To cos ^ I 

^(S,T)=0- YoCOS^ E2 TeCOsf ' 

V To COS f FeCOsf Ei J 

with the eigenenergies En ^ en + Un{n — l)/2 of the 
decoupled dot. The lowest energy E(^g x)=o = E2 + z 
then follows from the smallest root of the quartic equa- 
tion n±(-z^-2zC/-(re±ro)2cos2f) = {Eiz/2y. 
(ii) The {S,T) — 1/2 sector can be decom- 
posed into four subspaces with one or three 
electrons according to 5^ — ±1/2 and Cooper 
pair channel i> — e, o. The Hamiltonian is 

ii-(-) - ( ^1 cos f \ (e) 

^(S,T) = l/2 - ^^r^COsf E-i y ^^"^^^ "-{S.T) = 1I2 

operates in the subspacc spanned by {(ij^.j.|0), 6|(ij^.j,|0)} 
for S, = +1/2, and {d\^\S) ,b\d\^\G)} for ^ -1/2. 



FIG. 3: (Color online) (a) Same as Fig.[2]but for A = lOP 
within the zero-bandwidth limit for the leads (see text). 
Although the A — >■ 00 phase diagram is basically re- 
produced, for finite A, the {S,T) ~ 1/2 phase (green) 
exhibits a crossover from 0- to yr-junction behavior for 
f/ ~ A, as illustrated in panel (b), where the CPR is 
shown for t/T — —5 and several U ; the current is normal- 
ized to Jo, see Eq. @. Moreover, a phase with (5*, T) — 
at ((3 = and (S, T) = 1/2 ai ip — tv appears, where (con- 
trary to the "black" phase) <p = n corresponds to the 
lowest energy (7r'-behavior) , indicated in red [within the 
dashed ellipses in panel (a)]. 

Phase diagram for A ^ F. — Next we discuss the re- 
sulting phase diagram in the S'C/(4) limit {6 = 0). The 
result for A — >■ 00 is shown in Fig. [5] in the [/ — e 
plane. The phases are classified according to the three 
sectors defined above<^ The reported phases are specific 
for the SU{A) symmetry and are qualitatively different 
from the standard SU{2) case. We observe that the 
(ys-dependence of the A — > 00 ground-state energy im- 
plies 0-junction behavior for both S — and S = 1/2. 
While the magnetic S* = 1/2 sector often represents a 
7r-junction)^i^i^ in multi-level dots there is no direct con- 
nection between the spin and the sign of the Josephson 
couplingJ^ The 7r-phase found under perturbation the- 
ory [Eq. (|4]) for e < 0] is in fact restricted to the regime 
L/ > A, while for [/ < A, the S" = 1/2 state dis- 
plays a 0-phase. In the intermediate regime one should 
therefore observe a crossover between those two behav- 
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iors. Interestingly, there are parameter regions with a 
spin/pseudo-spin transition as (p varies. For instance, 
the "black" regions in Fig. [2] correspond to a mixed state 
with (5,r) = at ^ = and {S,T) = 1/2 at = tt, 
while for the "blue" region, the ground state is in the 
(S", T) = sector except at (/s = tt where it crosses to the 
(5, T) = 1 sector. 

We find that these phases are also observable at finite 
^ ^ r, where we have employed two complementary 
approaches. First, a full numerical solution is possible 
when approximating each electrode by a single site (zero- 
bandwidth limit), which can provide a satisfactory, albeit 
not quantitative, understanding of the phase diagram^ 
Second, one can go beyond the above A — >■ oo limit 
by including cotunneling processes in a systematic way. 
Both approaches give essentially the same results, and 
here we only show results from the single-site model. As 
can be observed in Fig. Ela), the overall features of the 
A — >■ oo phase diagram are reproduced for finite A, with 
somewhat shifted boundaries between the different re- 
gions. In particular, in the "green" (5* = 1/2) regime, 
this calculation captures the mentioned transition from 
a 0-junction at A L/ to a 7r-junction at A <C ?7, as 
illustrated in Fig. [3lb). Consequently, for finite A, the 
"black" phase may now have lowest energy at (ys = tt, im- 



plying the Tr'-phase^i^i^ indicated in "red" in Fig. [3ja) . 
Finally, for the junctions with U /T = 10.5 and 11 in 
Fig.EJb), the ground state is realized at phase difference 
< (/3 < TT, which implies that time-reversal symmetry is 
spontaneously broken here. 

To conclude, we have studied the Josephson current in 
SU{A) symmetric quantum dots, including the crossover 
to the standard SU{2) symmetric case. Contrary to 
normal-state transport, the supercurrent is very sensitive 
to the symmetry group, and should allow to observe clear 
signatures of the S'C/(4) state in ultra-clean CNT dots. 
In particular, the 7r-phase is largely suppressed, the CPR 
in the Kondo limit has a distinctly different shape and a 
smaller critical current, and the phase diagram turns out 
to be quite rich. In addition, following Ref. Ull, we ex- 
pect a strongly reduced thermal noise in the deep SU (4) 
Kondo regime since (in contrast to the SU{2) case) there 
are two channels with imperfect transmission. Future 
theoretical work is needed to give a quantitative under- 
standing of the crossover between the various regimes 
discussed above. 
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